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Abstract 

Density dependent relativistic Hartree-Fock theory has been extended to 
describe properties of exotic nuclei. The effects of Fock exchange terms and 
of vr- and p- meson contributions are discussed. These effects are found to 
be more important for neutron rich nuclei than for nuclei near the valley of 
stability. 
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Among the various theoretical models for the bulk properties of nuclei, the 
relativistic mean field theory (RMF) has been particularly successful [[l|, ^ |^]. 
Relativistic mean field theory can be regarded as an effective field theory which 
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achieves a description of nuclear properties throughout the entire periodic table 
with high precision. It is based on a relativistic Lagrangian which includes the 
nucleon, the omega-, the sigma-, and the rho-meson (which reproduces the asym- 
metry energy of nuclear matter) as well as non-linear self-interactions of scalar 
mesons. The inclusion of the non-linear self-interactions of the sigma field intro- 
duces density-dependent interactions in the model, which are crucial to reproduce 
quantitatively the bulk properties of nuclear matter and finite nuclei. Sharma et al. 
[01, Lalazissis et al. Maharana et al. [||], Hirata el al. [0] and Lalazissis et al. 
[Hi have systematically investigated the binding energies, deformations, rms radii, 
one (or two) nucleon separation energies, quadrupole moments and isotopic shifts 
of nuclei far from the stability line. A review of the relativistic mean field model 
has been given by Ring [^. 

Relativistic effects have been found to be important for theories starting from 
the realistic nucleon-nucleon interaction, too. The Dirac- Brueckner-Hartree- 
Fock(DBHF) approach based on the meson-exchange potential Bonn A success- 
fully reproduces nuclear matter saturation properties [[T0[|, while non-relativistic 
models failed. The DBHF approach has been applied to finite nuclei by Muether 
et al. [[TTt|. Due to the complexity of the DBHF in finite nuclei, only light spherical 
nuclei have been investigated. 

Brockmann and Toki have suggested a method that allows to link the binding 
energies and radii of finite nuclei to the nucleon-nucleon interaction via a local 
density dependence approximation to the DBHF [[T2|], the "relativistic density de- 
pendent Hartree" model(RDH). In this approach, density dependent coupling con- 
stants of the meson-nucleon interactions are derived from DBHF calculations in 
nuclear matter. Those coupling constants are used in a relativistic Hartree calcula- 
tion performed for finite nuclei. The model therefore is parameter free in the sense 
that no parameters are adjusted to data of the nuclear many-body problem. It does, 
however, contain some implicit assumptions. One assumes that it is sufficient to 
include density dependent modifications of the coupling constants of the sigma- 
and the omega-meson only. This assumption may be justified by the very success 
of the relativistic mean field model which shows that the sigma- and the omega- 
meson dominate the saturation properties of nuclei, as is obvious in closed-shell 
nuclei and in the Hartree approximation. Moreover, the Dirac Brueckner Hartree- 
Fock calculation provides a scalar and a vector contribution to the self-energy of 
the nucleon for each density. These are two constraints which can be used to de- 
termine two parameters of the effective Lagrangian. It is then natural to adjust the 
coupling constants of the two most relevant mesons to those numbers. 

The RDH approach with an effective sigma, omega, and rho meson gives satis- 
factory descriptions of ground state properties of a few stable nuclei [|T^, although 
the fit quality cannot match the one achieved in the RMF theories. Investigations 
of this type are very important for our understanding of properties of finite nuclei 
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in a parameter free way. The model so far has been applied only to a relatively 
small number of nuclei. Therefore it is of great interest to extend such density 
dependent relativistic model to deformed and exotic nuclei. Recently, the RDH 
has been applied to deformed nuclei and the investigations show large deviations 
from the experimental data for the quadrupole moments of rare-earth nuclei [[T^] . 

In the RDH approach, Pauli exchange effects have been neglected. Fritz, 
Miither and Machleidt [ [T4| ] pointed out, however, that the nuclear radii are re- 
duced once the Pauli exchange effects are taken into account. On the other hand, 
after incorporating both pion and rho-exchange (including in particular the tensor 
part of the rho-nucleon coupling) into the relativistic density dependent Hartree- 
Fock approach, larger radii are found [|T5|]. The findings of [[1^] suggest that the 
effects of Fock terms be investigated more thoroughly. In the case of spin-orbit 
splitting, the inclusion of the Fock terms appears to be essential. In "^^Ca, rela- 
tivistic mean field theory gives a splitting of the ^5/2 and the ^3/2 levels of 9.32 
MeV, whereas the relativistic density dependent Hartree-Fock theory gives 3.19 
MeV which comes closer to the experimental value of 4.3 MeV [[T5|]. 

Another place where the influence of Fock terms may be manifest is nuclei 
with a large neutron excess. The masses and the density distributions of neutron- 
rich nuclei are important to predict the abundance of elements [[T^. Systematic 
experimental investigations of exotic nuclei have become possible because of the 
recent advent of radioactive beam facilities [[TTp . 

In a study of neutron-rich nuclei, one has to include both the pion and the rho- 
meson. This requires to consider the effects of short-range correlations. Recently, 
Marcos et al. [[TSp have developed a nuclear matter model of the two-body cor- 
relation function. It provides a physical way to deal with the zero-range piece of 
TT and p interaction in the medium. In our investigation, this effect of short-range 
correlations is incorporated by removing contact interactions in tt and p exchanges 
[[T^p. In addition, Marcos et al. show that the effects of short range correlations 
can be simulated by modifiying the meson masses in the medium, at least in the 
Hartree approximation and for vanishing momentum transfer. 

In this letter, we apply the relativistic density dependent Hartree (RDH) and 
Hartree-Fock (RDHF) to stable as well as neutron rich nuclei. We extend the 
calculation to include the pion and the rho-meson using the following interaction 
Lagrangian: 

-Ci = 'ip{ga(T - Qu^l^ujf,- gpi^p,,- t)iIj (1) 

-V'e7^A^^(l + t^)^Ij - I^^^^^d^Ti ■ Tip 
where Tj indicate the isospin Pauli matrices. The effective strengths of couplings 
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Figure 1 : The binding energies of the Ca isotopes are shown. The open circles 
represent the experimental data and the full line with the filled circles show a 
fit using a liquid drop model I^^J. The RDH calculation including the rho meson ( 
dashed line with filled triangles ), the RDHF model with sigma and omega mesons 
only ( full line with open squares ), and the RDHF model including the isovector 
mesons ( full line with filled squares ) are compared. An RDH calculation without 
isovector mesons ( dashed line with open triangles ) is given as a reference point 
for the RDHF calculations. 



between the nucleon and mesons are denoted by the coupling constants gi or fi 
(i = a, cu, p, tt), respectively. The formalisms of the RDH and RDHF have been 
presented in detail in Refs. [[T^, |2l], [T^, [T^, p2[ ]. 

As Lenske pointed out rearrangement terms have to be included once density 
dependent coupling constants are used [pO|]. The inclusion of the rearrangement 
terms in the RDH slightly changes the binding energy, but improves the root- 
mean- square radii and density distributions. For simplicity, the rearrangement 
terms have not been included in this letter, since we are interested in the contribu- 
tions of the Fock term and isovector mesons, pion and rho in neutron-rich nuclei 
qualitatively. 

For the effective interaction, the coupling constants of the sigma- and omega- 
mesons are adjusted to reproduce the scalar and vector potentials and Uq in 
nuclear matter of the DBHF results using the Bonn A potential [|T^. The density 
dependence of the coupling constants is then adjusted to reproduce the nucleon 
self-energies at each density in the cases of RMF or RHF with or without isovector 
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Figure 2: The proton radii of the Ca isotopes are shown for some of the models of 
Fig.l. The calculations are denoted by the same symbols as in Fig.l. In addition 
the open circles show the experimental proton radii of^OCa and ^OCa. 

mesons. The nucleon and u meson masses are chosen to be the same as in the 
DBHF calculation, where M = 938.9 MeV, = 550 MeV and = 782.6 
MeV. For the ttNN vertex pseudo-vector coupling and for pNN vector and tensor 
couplings are used. The masses and coupling constants of isovector mesons are 

fixed to be m^ = 138 MeV, m^ = 770 MeV, g = 0.08, |f = 0.55 and ^ = 3.7 . 
In the RDHF calculation with isovector mesons, the sigma meson mass has been 
adjusted to m^r = 450MeV ( after removal of the contact terms of the pion and the 
rho meson ). The RDHF calculation without isovector mesons has been performed 
with rria = 550MeV. 

Our calculations are based on the one-boson exchange potential Bonn A 
which is used as an input for a Dirac Brueckner Hartree-Fock calculation. We will 
study the properties of exotic nuclei using three different approaches: 

1. Relativistic density dependent Hartree theory (RDH) with the a, to, and p 
mesons. 

2. Relativistic density dependent Hartree-Fock theory (RDHF) with a and to 
mesons only. 

3. Relativistic density dependent Hartree-Fock theory with a, u, n and p mesons 
(RDHF-i-TT-i-p). 
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Figure 3: The neutron radii of the Ca isotopes are shown for some of the models 
of Fig. 1. The calculations are denoted by the same symbols as in Fig.l. 

In Fig.l, we show the binding energies of the Ca isotopes obtained in our 
model. Experimentally, the binding energies of the Ca isotopes are known only 
till A=56 (open circles ) p3|]. We use the finite range liquid drop model by MoUer, 
Nix and Kratz [ p4| ] to extrapolate the binding energies of the exotic nuclei till 
A=70. The RDH approach allows to obtain a fair agreement with the experimental 

„2 

data, if one incorporates the rho meson with a coupling constant = 0.978) that 
has been chosen to reproduce the asymmetry energy. 

In both the relativistic density dependent Hartree- and Hartree-Fock- approxi- 
mation, the omega-nucleon and the sigma-nucleon coupling constants are adjusted 
to reproduce the nucleon self energy in nuclear matter. Therefore it is not too 
surprising that both methods reproduce the binding energies of N=Z nuclei with 
comparable accuracy. Indeed, Fritz et al. find in ^°Ca a binding energy per par- 
ticle of E/A = -8.21MeV in the RDH approach, while the RDHF calculation 
yields E/A = —7.76MeV which is in fair agreement with the experimental value 
of -8.50 MeV. The inclusion of the pion did not change the radius of "^^Ca, but 
slightly increased the binding energy [|T4l]. 

With increasing neutron excess, however, the inclusion of the Fock terms leads 
to a reduction of the binding energy. The inclusion of the pion and the rho leads 
to further decreases of the binding energies for the exotic nuclei. The comparison 
with the experimental data shows that all approaches: RDH with proper rho cou- 
pling and RDHF with and without isovector mesons could reproduce the slope of 



6 



1,E+00 

1,E-01 

1,E-a2 

i 1,E-03 
o 

1,E-04 
1,E-05 
1,E-06 

2 4 6 8 10 12 14 

r[fm] 



Figure 4: Neutron densities for "^"Ca: thin solid line: RDHF with uj and a only; 
thick solid line RDHF including isovector mesons; dashed line (hidden beneath 
the thick solid line) : RDHwith p. 

the curve of the experimental binding energies. At the RDHF level, the coupling 
constant of the rho-meson is not adjusted to the asymmetry energy, but taken from 
the underlying two-body interaction. 

The proton radii are displayed in Fig. 2a. In comparison with a RDH calcu- 
lation ( including the rho meson), the Fock term reduces the proton radii for all 
nuclei of the isotope chain. But the inclusion of the pion and in particular the rho- 
meson in the RDHF calculation increases the proton radii. For the neutron radii 
( Fig.2b ), a reduction due to the Fock term is found for mass numbers less than 
A=66, while for the extremely neutron rich nuclei, one even finds an increase of 
the neutron radii. The inclusion of the pion and the rho leads to an overall increase 
of the neutron radii. 

The various effects can be very clearly seen in the density distributions. We 
show here as examples the neutron density distributions for "^^Ca (Fig. 3a) and 
™Ca (Fig. 3b). We see that in the case of ^°Ca neither the Fock terms nor the 
inclusion of isovector mesons change the neutron density distribution consider- 
ably, so that one may argue in this case that the RDH description is sufficiently 
good. On the other hand we see strong effects of the Fock terms as well as the 
isovector mesons in the case of ™Ca. This indicates that especially in the case of 
very neutron rich nuclei, where e.g. a neutron halo is expected, one has to take 
into account all additional effects. To conclude, the density dependent relativistic 
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Figure 5: Neutron densities for ™Ca: thick solid line: RDHF including isovector 
mesons; thin solid line: RDHF with oj and a only; dashed line: RDH with p. 

Hartree-Fock theory has been employed to investigate properties of neutron rich 
nuclei. For the first time the effect of Fock exchange terms and isovector mesons 
on these properties has been studied. It is found that these terms give large contri- 
butions to properties of exotic nuclei. 
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